On the Detection of Magnetic Helicity 
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Magnetic fields in various astrophysical settings may be helical and, in the cosmological context, 
may provide a measure of primordial CP violation during baryogenesis. Yet it is difficult, even 
in principle, to devise a scheme by which magnetic helicity may be detected, except in some very 
special systems. We propose that charged cosmic rays originating from known sources may be useful 
for this purpose. We show that the correlator of the arrival momenta of the cosmic rays is sensitive 
to the helicity of an intervening magnetic field. If the sources themselves are not known, the method 
may still be useful provided we have some knowledge of their spatial distribution. 
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Magnetic fields pervade all astrophysical objects 0, Q 
and there are good theoretical reasons to believe that a 
weak magnetic field is present throughout the universe. 
In astrophysical systems the magnetic field is often he- 
lical which means that the field lines are twisted (like 
corkscrews), or that closed magnetic lines are linked. 
Mathematically, the average helicity density in a volume 
V is defined as: 
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H =— d^x A-B. 
V Jv 

In cosmology, a number of scenarios predict the creation 
of a primordial field with non-zero helicity. In the sce- 
nario discussed in Refs. jEH^IEISl; a magnetic field is pro- 
duced at the electroweak phase transition. The helicity 
of the magnetic field is related to the cosmological baryon 
asymmetry arising from CP violation in the fundamental 
particle physics theory, and the sign of the helicity is pre- 
dicted to be left-handed . There are also several other 
scenarios for the generation of primordial helical mag- 
netic fields that do not depend on the dynamics through 
a phase transition 0,lllallSinilIlll3- 

The helicity of magnetic fields in astrophysical jets can 
be deduced from the polarization of synchrotron radia- 
tion 14, JjJ. In such situations, the velocity of electrons 
in the jets is known and this additional information is 
crucial to the determination of helicity. In other situa- 
tions, it is much harder to find the helicity. For example, 
Faraday rotation only provides an estimate of the line of 
sight component of the magnetic field. Even by observing 
the Faraday rotation from different sources, the informa- 
tion is insufficient to estimate the helicity 0, IT^ . 
An estimate of the helicity necessarily requires sensitiv- 
ity to all components of the magnetic field and hence 
it is a challenging theoretical problem to devise means 
by which it may be measured. In Ref. ,2^ the 
imprint of cosmological magnetic helicity in parity-odd 
cross correlations of the cosmic microwave background 
(CMB) fluctuations was investigated, while in Ref. 
it was shown that helicity would introduce circular polar- 
ization of induced relic gravitational waves. Both these 



potential signatures of helicity are limited to cosmological 
magnetic fields since they rely on properties of the cos- 
mic microwave background or on the cosmic gravitational 
wave background. Further, the signals are small for sev- 
eral reasons: the cosmic magnetic field is constrained to 
be weaker than ~ 1 nG, the CMB is polarized only at 
the 10% level, the tensor modes that enter parity-odd 
correlations are tiny, and the gravitational waves are ex- 
tremely weak. 

In the present paper, we show that correlators of the 
arrival momenta of charged cosmic rays from known 
sources carry information about the helicity of the mag- 
netic field through which the charges propagate. The 
scheme has the advantage that it is not restricted to 
cosmological magnetic fields, and it utilizes cosmic rays 
which are abundant and well-studied . The difficulty 
with our scheme is that we do not normally know the 
source from which an observed cosmic ray emanated. 
However, the scheme may be extended to situations 
where we have some knowledge of the distribution of 
sources e.g. if the sources are located within a certain 
region of the disk of the galaxy. We have not yet ex- 
plored this possibility in detail. For the present paper, 
we focus on establishing an observable that is sensitive 
to magnetic helicity. Further work is needed to decide if 
the observable that we propose is practically useful. 

Consider a situation where there are two known sources 
(A and B) that are emitting charged particles that arrive 
on Earth. The particles would propagate along straight 
lines from the sources to the Earth if there were no mag- 
netic field. However, the trajectories get bent by the 
weak magnetic field. We work to lowest order in the 
magnetic field strength and consider the momenta of the 
particles as being perturbed due to the magnetic field: 



A = 



OA 



PA 



B = 



OB 



Pb 



where the in the subscript denotes an unperturbed mo- 
mentum, and pa.b are the momentum perturbations. 
The unperturbed momenta are directed along the lines of 
sight to the sources and the magnitudes are completely 
determined by the energies of the charged particles. 
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FIG. 1: Two sources A and B emit charged particles that are 
observed on Earth at O. If there was no ambient magnetic 
field, the particles would follow straight trajectories (dashed 
lines). In the presence of a weak magnetic field, the trajecto- 
ries get bent (solid curves). 



We are interested in the observable 

{PA{tf)Ph{t'f)) - PoaPoB + {PA{tf)pkt'f)) (1) 

where = 1,2,3 are spatial indices and tf and t'^ de- 
note arrival times from the two sources. The ensemble 
average refers to an average over many realizations of the 
magnetic field for the same locations of the two sources. 
We will discuss ways in which an ensemble average can 
be implemented toward the end of the paper. In writ- 
ing Eq. ^ we have implicitly considered particles for 
which the energies are fixed. Otherwise we would also 
have to average the unperturbed momenta since these 
depend on the energies of the particles. We have also 
taken (pa,b) = which holds for a stochastic magnetic 
field with zero mean. 

The first term in Eq. contains the unperturbed 
momenta. To evaluate this contribution, it is essential 
that we have some knowledge of the locations of the two 
sources. We begin by assuming that we know the loca- 
tions exactly and later comment on the case where the 
sources are distributed in a plane. We now evaluate the 
correlator for the momentum perturbations 

c^^' ^c^'O^A,^B)^{Mtf)pi{t'j)). 

where and are the vectors from O to sources A 
and B respectively. We find that certain components of 
C" are sensitive only to the helicity of the intervening 
magnetic field and vanish if the helicity is zero. 

To proceed, we use the Lorentz force law to find p to 
linear order in the magnetic field B, 



q 



^ dt' vo X B(x(i')), 



(2) 



where we have temporarily suppressed the subscripts 
specifying the source for convenience. The unperturbed 
velocity vq may be related to the unperturbed momen- 
tum using Vq = Fq/E, where E is the energy. The initial 



momentum perturbation pi = p(ti) need not vanish and 
also depends on the intervening magnetic field. To de- 
termine Pi we integrate P = Pq + p(i) over proper time 



x(0 



E 



E 



dt' dt" Vo X B(x(t")), 



(3) 



where Xi = x(ii) and E is the energy of the particle. 

We know that the charged particle arrives at the detec- 
tor aXt — tf. Taking the origin of the coordinate system 
to be at the position of the detector, we find 

p, = -| dt' ^' dt" vo X B(x(t")) (4) 

where T = t f — ti and we have used Po = —Ex{ti)/T. 
Going back to Eq. jSJ, we can write 



p(tf) = J[vo X B(x)] 



(5) 



where the action of the operator 2 on a function F(t) is 
defined by 

I[F] = J\t' J* dt"F{t") + q J'' dt'F{t') 



T 



^ / dt'[l-T5{t' -tf)] / dt"F{t"){&) 



Therefore 



^^r,ue,rk'vlj,vls{B\^A{t))B^'{^B{t'))) 



(7) 

where, for generality, we consider different species of 
charged particles arriving from sources A and B with 
different energies Ea, Eb and travel times Ta, Tb- 

The auto-correlator of an isotropic, stochastic, time- 
independent magnetic field can be written as 



(B,(x + r)B,(x)) - Mn{t) 



+ MH{r)eijiri, 



(8) 



where MN{r), ML{r), and Af^f(r) are the correlation 
functions for the "Normal" , "Longitudinal" , and "He- 
lical" parts of the magnetic field. We assume that space- 
time curvature can be neglected on the length scales of 
interest and use a Minkowski metric, so Bi = — S*. All 
correlation functions depend only on r = |r|, reflecting 
the statistical isotropy of the field. The divergence-less 
condition requires 

The ensemble averaging in Eq. © is over all locations x 
but for fixed r. 
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The magnetic field two-point correlation function is 
often given in Fourier space, so it is useful to express 
Mjv(r), Mi(r), and Mnir) in terms of a magnetic field 
wave number k space power spectrum defined by: 



kikj, and the unit 



where the projector Pij(k) = 6 
vector ki = ki/k. F]\[{k) and Fnik) are the symmetric 
and helical parts of the magnetic field power spectrum, 
related to the average energy density and helicity of the 
magnetic field. The functions Fat (A) and Fnik) can be 
related to the correlation functions MN{r), ML(r) and 
Mnir) as in Ref. 

The correlator C" , Eq. may be decomposed into 
normal, longitudinal, and helical parts. 



C" + cf + c 



H 



The remaining calculation involves inserting Eq. (jS)) into 
((TJ and simplifying. Let us define 



n = voA X vos • 



(9) 



A straightforward computation gives the correlator in- 
duced by the normal component of the magnetic field 
power spectrum. 



MN{r) (5" voA • vos - Vq^v, 



OB 



-[voA X rJ'[voB X rj 



(10) 



where r{t,t') — xyi(t) — xb(<') and the unit vector f = 
r/r. The longitudinal piece of correlator is 



MLir)[voA X f]'[voB X f]' 



(11) 



Similarly for the helical component we get: 



Mnir) { e"''[(voA • r)i;[,B + «(,^(vob • r)] 



Cj/ — ^A^B 



The helical part of correlator, C|| , vanishes for i = i', 
and the trace of the momentum perturbation correlator 
contains contributions only from the normal and longi- 
tudinal parts of the magnetic field spectrum: 



C 



TR 



IaIb 



2MAr(r)voA • VoB 



-(Mi(r) - MN{r))[{^roA ■ ^ob) - (vqa • f)(voi3 • f)] 



Let us take our coordinate system so that the triangle 
ABO lies in the xy— plane (see Fig. ^ and n is in the 
z— direction. We find that all components of the helical 
correlator vanish except for: 



H ~ 

r^zi 

— 



-2(voB •Mh)(voa X %f 
-2(voA •MH)(voi3 X z)* 



where 



(12) 
(13) 



(14) 



After doing the integrations, any dependence of M/f on 
and i!j can be traded for a dependence on the posi- 
tion of the particles at the final time using: X/i(</) = 
xs(i'y) = 0. A neat combination of the two components 
in Eqs. (O and lO is 



C 



Hi 



M'^lvoA X vosi 



(15) 



where we have used Mh ■ z = 0. 

The normal component C"j^ vanishes if i or i' (but not 
both) are in the z— direction. The trace of momentum 
correlator has the contribution only from Mn and Ml, 
and in a;?;— plane = i'^ ~ ^' v) depends only 

on the normal component MN{r), 



= lAlB[MNir)] VoA • VOB 



The longitudinal component C" vanishes if i ^ i' and 
has only one non-zero component when i = i' is along 
the z direction. 



Cf 



IaIb 



AfL(r){voA • voB - (vqa ■ f)(voyi • f)} 



On the other hand, C)} is non-zero only if one (and only 
one) of i, i' is along the z— direction. This can be under- 
stood on physical grounds as follows. If the magnetic field 
is not helical, a charged particle is as likely to be deflected 
in the +z direction as it is to be in the —z direction by 
the stochastic magnetic field. By symmetry, the compo- 
nents C*^ must then vanish. However, a helical magnetic 
field breaks the symmetry and these components become 
non-zero. So we see that only the helical contribution 
enters the xz, yz, zx and zy components of C" , and fur- 
ther, only the non-helical contributions enter the other 
components. Therefore the normal and helical pieces of 
the correlator do not mix. 

The other components of the correlator {e.g. C'^y) can 
be used to find the normal and longitudinal correlation 
functions, and AI^. Correlations of the rotation mea- 
sure due to Faraday rotation of polarized light from dif- 
ferent sources can also be used as an independent method 
to determine Mjv and Ml [l^ - 

The above analysis relies on an average over an ensem- 
ble of random magnetic field realizations for fixed source 
and detector positions (triangle ABO in Fig. In our 
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universe, however, only one realization of the magnetic 
field is available and so we need to discuss a practical 
scheme for doing the ensemble average and thus estimat- 
ing the correlator, C" . The ensemble average would 
necessarily involve averaging over many pairs of sources 
that we denote by {A,B)a, a — 1,2,...,N from which 
we assume that cosmic rays are observed with average 
momenta (P,P')a- Going back to the magnetic field 
correlation function, Eq. we see that the ensemble 
average is over all x for fixed point separation r. Since 
C" depends on integrals of the magnetic field along the 
line of sight, this suggests that we find pairs of sources 
such that r(t,t') is the same function of t and t' for all 
of them. Together with the constraint that the posi- 
tion of the detector is fixed at O, such an ensemble has 
only one element and is not useful. Instead of holding 
r{t,t') fixed, it is more useful to choose pairs with a less 
restrictive condition. There are many such possible con- 
ditions, each with its own advantage. As an example, 
one such condition is that the source separation vector, 
A = Xs — Xyi, be held fixed for all pairs in the ensemble. 
The corresponding observable is 

(A) = - 5] (16) 

This observable quantity is an estimator for an average 
of the total correlator C" discussed above, where the 
average is over sources with fixed separation vector A. 
As we have noted before, the helical part, and only the 
helical part, of the magnetic field enters certain compo- 
nents of the total correlator (see Eqs. (|12(l and l|13() '). So 
to get an estimator for the helical components of the cor- 
relator we should only look at those components of the 
estimator that involve momentum in a direction in 
the source-observer plane and the second momentum in 
a direction perpendicular to this plane. Since different 
pairs of sources (with the same A) may lie in different 
planes, an estimator for the helical part of the correlator 
is 

1 ^ 

a=l 



where Hq, is the normal to the plane of the source pair 
labeled by a and the observer as defined in Eq. and 
is some chosen unit vector within the plane. 

To extract the magnetic helicity from a measurement of 
Ph{^), wc must evaluate the corresponding theoretical 
quantity, which is given by 

CH(A) = i / d3x^Cff(X^,XA + A) 

y Jv 

The above integral is quite involved but can be done nu- 
merically for different choices of M// (r) . 

To summarize, a non- vanishing observed value of 
■Pij(A) will give a measure of C_y(A) and hence will 
lead to the magnetic helicity, M//(r). 

At present we do not have any known sources of 
charged cosmic rays. However, it is likely that a large 
fraction of cosmic rays that we see arise in the galactic 
disk (or sources confined in the cosmic large-scale struc- 
ture 123 )• It may be possible to usefully extend the en- 
semble average to include pairs of locations in the galac- 
tic disk. Such an averaging could still yield information 
about the helicity of the galactic magnetic field. We plan 
to consider this extension of our result in future work, 
together with other observational issues. 
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